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Context. A phase of strong interacting matter with deconfined quarks is expected in the core of a massive neutron star. If this decon- 
finement phase transition is of the first order, as suggested by many models inspired by quantum chromodynamics, then it will be 
QQ ' triggered by the nucleation of a critical size drop of the (stable) quark phase in the metastable hadronic phase. Within these circum- 

stances it has been shown that cold (T = 0) pure hadronic compact stars above a threshold value of their gravitational mass (central 
pressure) are metastable with respect to the "decay" (conversion) to quark stars (i.e., compact stars made at least in part of quark 
matter). This stellar conversion process liberates a huge amount of energy (a few 10^^ erg), and it could be the energy source of some 
' of the long gamma ray bursts. 

• ' Aims. The main goal of the present work is to establish whether a newborn hadronic star (proto-hadronic star) could survive the early 

i-S^ ' stages of its evolution without "decaying" to a quark star. To this aim, we study the nucleation process of quark matter in hot (T i= 0) 

/J-stable hadronic matter, with and without trapped neutrinos, using a finite temperature equation of state (EOS) for hadronic and 
quark matter. 

' Mettiods. The finite-temperature EOS for the hadronic and for the quark phases were calculated using the nonlinear Walecka model 

, and the MIT bag model, respectively. The quantum nucleation rate was calculated making use of the Lifshitz & Kagan nucleation 

^ , theory. The thermal nucleation rate was calculated using the Langer nucleation theory. 

Results. We calculate and compare the nucleation rate and the nucleation time due to thermal and quantum nucleation mechanisms. 
We compute the crossover temperature above which thermal nucleation dominates the finite temperature quantum nucleation mech- 
anism. We next discuss the consequences of quark matter nucleation for the physics and the evolution of proto-neutron stars. We 
^ , introduce the new concept of limiting conversion temperature and critical mass M^,- for proto-hadronic stars, and we show that proto- 

1/^ ■ hadronic stars with a mass M < M„ could survive the early stages of their evolution without decaying to a quark star. We extend the 

' concept of maximum mass of a "neutron star" with respect to the classical one introduced by Oppenheimer & Volkoff to account for 

I the existence of two distinct families of compact stars (hadronic stars and quark stars) as predicted by the present scenario. 
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1 . Introduction order phase transitions are triggered by the nucleation of a crit- 
ical size drop of the new (stable) phase in a metastable mother 

According to quantum chromodynamics (QCD) a phase transi- phase 
tion from hadronic matter to a deconfined quark phase should 

occur at a density of a few times nuclear matter saturation den- in j-gcent years there has been a growing interest in study- 

sity po ~ 2.8 X lO'-* g/cm^ . Neutron star structure calculations, the nucleation process of quark matter in the core of mas- 

■ - - based on a wide variety of modern equatio ns of sta te (EO S) of ^ive neutron stars. The nucleation rate of quark matter in cold 

hadronic matter (Lattimer & Prakash|20(JlJ Bombaci|2007}, pre- (7^ ^ 0)_^stable nuclear matter has been calculated by lida & 

diet a maximum stellar central density (the one for the maxi- Sato (fT997] | wthin the Lifshitz-Kagan quantum nucleation the- 

mum mass star configuration) in the range of 4 - 8 times po- ory (L ifshitz & Kagan ,1912). The same authors (lida & Sato 

Consequendy, the core of neuti'on stars is one of the best candi- have later on studied quark matter nucleation in y6-stable 

dates in the universe where a phase of sti'ong interacting matter hyperonic matter within a relativistic extension of the Lifshitz- 

with deconfined up, down and strange quarks could be found. Kagan theory. Various significant astrophysical implications of 

In the region of high density (high baiyon chemical poten- quark matter nucleation in neutron stars have bee n explored in 

tial) and low temperature (which is the one relevant for neu- a series of subsequent papers (Berezhiani et al. |20521 |20031 

tron star physics) many QCD-inspired models suggest the de- Bombaci et al. 120041 Drago et al. l2004l L ugones & Bombaci 

confineme nt trans ition to be a first-order phase transition (Hsu 2005; Bombaci et al. '2007', Bombaci et al. '2008', Drago et al. 

& Schwetz 119981 Fodor & Katz 2004 ). As is well known, first- ^008 ; Bambi & Drago 2008). In particular, in these studies it has 

been shown that, above a threshold value of the central pressure, 

* corresponding author email: bombaci@df.unipi.it a pure hadronic compact star (HS) is metastable to the decay 
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(conversion) to a quark star (QS), i.e. to a hybrid neutron star 
or to a strange star (Bodmer 1971' Witten 19841, depending on 
the details of the EOS for quark matter used to model the phase 
transition. This stellar conversion process releases (Bombaci & 
Datta ,2000J an extraordinarily large amount of energy (a few 
10^^ erg) and it could be the energy source of some of the long 
gamma ray bursts (GRBs). 

The research reported in these papers has focused on the 
quark deconfinement phase transition in cold and neutrino-free 
neutron stars. In this case the formation of the first drop of QM 
could take place solely via a quantum nucleation process. 

A neutron star at birth (proto-neutron star) is very hot (T - 
10-30 MeV) with neutrinos being still trapped in the stellar 
interior (Burrow & Lattimer 1986; Prakash et al. |1997| Pons et 
al.[1999). Subsequent neutrino diffusion causes deleptonization 
and heats the stellar matter to an approximately uniform entropy 
per baryon § =1 - 2 (in units of the Boltzmann's constant feg). 
Depending on the stellar composition, during this stage neutrino 
escape can lead the more "massive" stellar configurations to the 
formation of a black hole (Bombaci. 1 1996t Prakash et al. |1997l l. 
However, if the mass of the star is sufficiently low, the star will 
remain stable and will cool to temperatures well below 1 MeV 
within a cooling time tc„oi ~ a few 10^ s, as the neutrinos con- 
tinue to carry energy away from the stellar material (Burrow & 
Lattimer [T9g6l Prakash et al. IT9971 Pons et al. [1999). Thus in 
a proto-neutron star, the quark deconfinement phase transition 
will likely be triggered by a thermal nucleation process. In fact, 
for sufficiently high temperatures, thermal nucleation is a much 
more efficient process than the quantum nucleation mechanism. 

Some of the earlier studies of quark matter nucleation (see 
e.g., Horvath et al. 1992] Horvath 1994; Olesen & Madsen 1994; 
Heiselberg 1995; Harko et al. 12004) have akeady dealt with ther- 
mal nucleation in hot and dense hadronic matter In these stud- 
ies, it has been found that the prompt formation of a critical size 
drop of quark matter via thermal activation is possible above a 
temperature of about 2-3 MeV. As a consequence, it was in- 
ferred that pure hadronic stars are converted to quark stars within 
the first seconds after their birth. However, these works (Horvath 
et al. 119921 Horvath 1994; Olesen & Madsen 1994; Heiselberg 
1 19951 1 reported an estimate of the thermal nucleation based on 
"typical" values for the thermodynamic properties characteriz- 
ing the central part of neutron stars. Quark matter nucleation 
during the early post-bounce stage of core collapse supernova 
has been lately examined by Sagert et al. (I20091 l and Mintz et al. 
(I20TOI 1. 

The main goal of the present paper is to establish whether 
a newborn hadronic star (proto-hadronic star) could survive the 
early stages of its evolution without "decaying" to a quark star. 
To this aim, we calculated the thermal nucleation rate of quark 
matter in hot (T + 0) y6-stable hadronic matter using a finite 
temperature EOS for hadronic and quark matter In addition, we 
calculated the quantum nucleation rate at finite temperature, and 
compared the thermal and quantum nucleation time at differ- 
ent temperatures and pressures characterizing the central condi- 
tions of metastable proto-hadronic compact stars. We computed 
the crossover temperature above which thermal nucleation dom- 
inates the finite temperature quantum nucleation mechanism. 
An exploratory study of this issue has been recently reported 
(Bombaci et al. 120091 ) for the case of neutrino-free hadronic mat- 
ter Here, we extend our analysis and our numerical investigation 
to the case of hadronic matter with trapped neutrinos, and we 
make a systematic study of the finite temperature nucleation pro- 
cess on the value of the Bag constant in the EOS for the quark 



phase. Finally, we discuss the relevant consequences of quark 
matter nucleation for the evolution of proto-neutron stars . 

2. Equation of state 

Within a fundamental approach, the EOS for neutron star mat- 
ter should be derived by solving numerically the equations of 
QCD on a space-time lattice. Lattice calculations at zero baryon 
chemical potential (zero baryon density) suggest that at high 
temperature and for physical values of the quark masses, the 
transition to quark gluon plasma is a "crossover" (Fodor & Katz 
2004 ; Karsch 2005 ) rather than a phase transition, which would 
be signaled by singularities in the thermodynamic observables. 
Unfortunately, present lattice QCD calculations at finite den- 
sity (baryon chemical potential) are plagued with the notorius 
"sign problem", which makes them unrealizable by all presently 
known lattice methods (see e.g. Lombardo (2007) and references 
therein). Thus, to explore the QCD phase diagram at low tem- 
perature T and high baryon chemical potential y., it is necessary 
to invoke some approximations in QCD or to apply some QCD 
effective model. In this region of the T-yU plane, many QCD in- 
spired models suggest the deconfinement transition to be a first- 
order phase transition (Hsu & Schwetz 1998. Fodor & Katz 
2004). In this domain of the QCD phase diagram, many pos- 
sible color superconducting phases of quark matter are expected 
(see e.g. Casalbuoni & Nardulli 2004; Alford et al. 2008 and 
references therein). The effects of these color superconducting 
phases on the properties of proto-neutron stars have been inves- 
tigated in (Aguilera et al. 2004; Sandin & Blaschke 2007; Gu et 
al. 120081 Lugones et al. 2009 ). 

Very recently, a new phase of QCD, named quarkyonic 
phase, has been predicted (McLerran & Pisarski 120071 Hidaka 
et al. '2008). This matter phase is characterized by chiral sym- 
metry and confinement. 

Here, we have adopted a more traditional view, assum- 
ing a single first-order phase transition between the confined 
(hadronic) and deconfined phase of dense matter, and we used 
rather common models for describing them. 

For the hadronic phase we used the nonlinear Walecka model 
(NLWM) (Walecka [T9741 Serot & Walecka [T986] l. This model 
is based on a relativistic Lagrangian of baryons interacting via 
the exchange of various mesons. We considered a version of the 
NLWM in which the baryon octet particles (n, p. A, Yr , 2", 2^, 
S", H") interact via the exchange of the scalar cr, the vector- 
isoscalar and the vector-isovector meson fields (see e.g. 
Prakash et al. 1 1 9971 Glendenning 2000. Menezes & Providencia 
l2003T l. 

The Lagrangian density of the model reads as 

— -^hadrons "t" -^leplons (1) 

where the hadronic contribution is 

-Chadrons — -^baryons "t" -^mesons (2) 

with 

-Cbaryons^ J] lj> D - M^] l// , (3) 

baryons 

where 

£>/i = id^ - gu,B(^^i - gpBts ■ Pii, (4) 

and M*g - Mb - gaBC. The quantity ts designates the isospin of 
baryon B. The mesonic contribution reads as 

-^mesons — -^tr "t" -^o) "t" -^p? (5) 
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with 



1 1 , 



(6) 



(7) 



For the lepton contribution we take 
leptons 



(9) 



where the sum is over electrons, muons and neutrinos for matter 
with trapped neutrinos. In uniform matter, we get for the baryon 

Fermi energy efB = gioBi^o + gpBhBPm + -^^^^ with the 

baryon effective mass M*g - M - ga-BO'- 

In the present work we used one of the parametrizations of 
the NLWM given by Glendenning & Moszkowski (2001) (here- 
after the GMl equation of state). The parameters of the GMl 
EOS are fitted to the saturation properties of symmetric nuclear 
matter: binding energy per nucleon B/A = -16.3 MeV, satura- 
tion density po = 0.153 fm"-', incompressibility K = 300 MeV, 
nucleon effective mass M* = 0.7, and to the nuclear symmetry 
energy at saturation density Osym - 32.5 MeV. 

Including hyperons involves new couplings, i.e., the 
hyperon-nucleon couplings: g^^B = Xo-b go-, gcoB = 
XojB ga>, gpB - XpB gp- For nucleons we take Xo-b, x^b, XpB - 1 
and for hyperons we will consider the couplings proposed by 
Glendenning and Moszkowski ( 20011) . They consider the bind- 
ing energy of the A in nuclear matter, B/^, 



(^) 



-28 MeV 



Xc^ goj (^0 - Xcr ga-Cr 



(10) 



to establish a relation between x^- and x^. Moreover, known neu- 
tron star masses restrict x^r to the range 0.6-0.8. We take Xp - x^r 
and consider x^r - 0.6. 

For the quark phase, we adopted a phenomenological EOS 
(Farhi & Jaffe 1984) that is based on the MIT bag model for 
hadrons and considered different possible values for the bag con- 
stant B. The remaining model parameters have been fixed to: 
nil, - = 0, - 150 MeV for the masses of the up, down, 
and strange quark, respectively, and as - Q for the QCD struc- 
ture constant. 

The two models for the EOS have been generalized to the 
case of finite temperature (see e.g. Menezes & Providencia 
l2003T l. 

3. Phase equilibrium 

For a first-order phase transition, the conditions for phase equi- 
librium are given by the Gibbs' phase rule 



IUH(T,Po)^fiQ(T,Po) 
where 



Ph^Pq = Po 



sq + Pq- sqT 
ng 



(11) 
(12) 

(13) 
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Fig. 1. (Color online.) Gibbs energy per baryon, in neutrino- 
free matter, as a function of pressure for the hadronic and quark 
phases at different temperatures. Lines with the steeper slope re- 
fer to the hadronic phase. Full dots indicate the transition pres- 
sure Po for each temperature. GMl EOS for the hadronic phase 
and MIT bag model EOS with B = 85 MeV/fm^ for the quark 
phase. 



are the Gibbs' energies per baryon (average chemical potentials) 
for the hadron and quark phase respectively, and su {bq), Ph 
(Pq), Sfj (sq) and «// (n^) respectively denote the total (i.e., in- 
cluding leptonic contributions) energy density, total pressure, to- 
tal entropy density, and baryon number density for the hadron 
(quark) phase. Above the "transition point" (Pq), the hadronic 
phase is metastable, and the stable quark phase will appear as a 
result of a nucleation process. 

Small localized fluctuations in the state variables of the 
metastable hadronic phase will give rise to virtual drops of the 
stable quark phase. These fluctuations are characterized by a 
time scale v^' ~ 10"^^ s. It is set by the strong interactions 
(responsible for the deconfinement phase transition), and it is 
many orders of magnitude shorter than the typical time scale for 
the weak interactions. Quark flavor must therefore be conserved 
during the deconfinement transition. We refer to this form of de- 
confined matter, in which the flavor content is equal to that of the 
y6-stable hadronic system at the same pressure and temperature, 
as the Q*-phase. Soon after a critical size drop of quark matter 
is formed, the weak interactions will have enough time to act, 
changing the quark flavor fraction of the deconfined droplet to 
lower its energy, and a droplet of y6-stable quark matter is formed 
(hereafter the Q-phase). 

This first seed of quark matter will trigger the conversion 
(OUntoT987; Heiselberg et al. 1991; Bombaci & Datta|2000l 
Drago et al. 2010) of the pure hadronic star to a hybrid star or to 
a strange star Thus, pure hadronic stars with values of the central 
pressure higher than Pq are metastable to the decay (conversion) 
to hybrid stars or to strange stars (Berezhiani et al. 120021 120031 
Bombaci et al. I2004t Drago et al. 120041 Lugones & Bombaci 
2005; Bombaci et al. 120071 ). The mean lifetime of the metastable 
stellar configuration is related to the time needed to nucleate the 
first drop of quark matter in the stellar center, and it depends 
dramatically on the value of the stellar central pressure. 

In Figs. 1 and 2, we plot the Gibbs' energies per baryon for 
the hadron-phase and the Q*-phase in neutrino-free matter, at 
different temperatures (T = 0, 10, 20, 30 MeV) and for two dif- 
ferent values of the bag constant B = 85 MeV/fm^ (Fig. 1) and B 



4 



Bombaci et al.: Effects of quark matter nucleation on the evolution of proto-neutron stars 



= 100 MeV/fm^ (Fig. 2). Lines with the steeper slope refer to the 
hadron phase. As we see, the transition pressure Pq (indicated by 
a full dot) decreases when the hadronic matter temperature is in- 
creased. 

The phase equilibrium curve PoiT) for the hadron-quark 
phase transition (in the case B = 85 MeV/fm^) is shown in Fig. [3] 
for neutrino-free matter and matter with trapped neutrinos . The 
region of the Pq-T plane above each curve represents the decon- 
fined Q*-phase. As expected (Prakash et al. 119971 Lugones & 
Benvenuto'1998', Vidana et al. 2()05', Lugones et a. 12009) neu- 
trino trapping in y6-stable hadronic matter inhibits the quark de- 
confinement phase transition, thus the global effect of neutrino- 
trapping is to produce a shift of the phase equilibrium curve to- 
ward higher values of the pressure in the Pq-T plane. 

As is well known, for a first-order phase transition the deriva- 
tive dPo/dT is related to the specific (i.e. per baryon) latent heat 
Q of the phase transition by the Clapeyron-Clausius equation 

dPo _ tiHtiQ- Q 
dT hq' — nn T 

Q^Wq.-Wh = T{Sq.-Sh) 



(14) 



(15) 



where Wh (Wq-) and S h (§ q-) denote the enthalpy per baryon 
and entropy per baryon for the hadron (quark) phase, respec- 
tively. The specific latent heat Q and the phase numbers den- 
sities riH and «q. at phase equilibrium are reported in Tables 
1 and 2 for the case of neutrino-free matter and for two dif- 
ferent values of the bag constant B = 85 MeV/fm^ (Table 1) 
and B = 100 MeV/fm^ (Table 2). As expected for a first-order 
phase transition, one has a discontinuity jump in the phase num- 
ber densities: in our particular case nQ-{T, Pq) > tinij, Pq). This 
result, together with the positive value of Q (i.e. the deconfine- 
ment phase transition absorbs heat), tells us (see Eq. (fT4l l) that 
the phase transition temperature (ITTl decreases with pressure (as 
in the melting of ice). 

To explore the effect of neutrino trapping on the phase equi- 
librium properties of the system, in Table 3 we report the quan- 
tities Q, hh, rtQ- and Pq at different temperatures, for a bag con- 
stant B = 85 MeV/fm-' (compare with the results in Table 1). As 
we see, the phase number densities hh and ng. at phase equilib- 
rium are shifted to higher values, and the specific latent heat Q 
is increased with respect to the neutrino-free matter case. 

We also made a systematic calculations of the phase equi- 
librium properties varying the value of the bag constant in the 
EOS of the Q*-phase. These results are summarized in Fig. 4 
where we plot the phase equilibrium pressure Pq between the 
hadron phase and the Q* phase at zero temperature as a function 
of B for neutiino-free matter and neutrino-trapped matter. As 
expected, a higher value of B makes the deconfinemnt transition 
more difficult, increasing the value of the transition pressure. 

4. Quantum and thermal nucleation rates 

Above the phase equilibrium pressure Pq the jS-stable hadron 
phase is metastable, and the formation of the stable (with respect 
to the strong interactions) Q*-phase will occur via a nucleation 
process. The main effect of finite temperature on the quantum 
nucleation mechanism of quark matter is to modify the energy 
barrier separating the quark phase from the metastable hadronic 
phase. This energy barrier, which represents the difference in the 
free energy of the system with and without a Q*-matter droplet, 
can be written as 



> 

s 




50 55 
P [MeV fm"" 

Fig. 2. S ame as Fig . 1 , but with a bag constant B = 1 00 Me V/fm^ . 
(Color online) 

Table 1. The specific latent heat Q and the phase number 
densities hh and «q. at phase equilibrium. GMl EOS for the 
hadronic phase, MIT bag model for the quark phase with B = 
85 MeV/fm^. Results for neutiino-free matter. 



Table 2 



T 


Q 


JlQ. 


'% 


Pa 


MeV 


MeV 


fm-3 


fm-^ 


MeV/fm^ 





0.00 


0.453 


0.366 


39.95 


5 


0.56 


0.451 


0.364 


39.74 


10 


2.40 


0.447 


0.358 


38.58 


15 


5.71 


0.439 


0.348 


36.55 


20 


10.60 


0.428 


0.334 


33.77 


25 


17.17 


0.414 


0.316 


30.36 


30 


25.44 


0.398 


0.294 


26.53 


2. Same as Table 1, but with E 


= lOOMeV/fml 


T 


Q 


llQ. 


'% 


Po 


MeV 


MeV 


fm-3 


fm-^^ 


MeV/fm^ 





0.00 


0.557 


0.447 


59.53 


5 


0.44 


0.556 


0.445 


59.22 


10 


1.98 


0.552 


0.441 


58.18 


15 


4.91 


0.544 


0.432 


56.39 


20 


9.35 


0.533 


0.421 


53.81 


25 


15.49 


0.519 


0.401 


50.47 


30 


23.51 


0.503 


0.386 


46.44 


3. Same as Table 1, but with trapped neutrinos. 


T 


Q 


llQ. 




Po 


MeV 


MeV 


fm-3 


fm-3 


MeV/fm^ 





0.00 


0.603 


0.516 


113.77 


5 


0.65 


0.601 


0.514 


113.11 


10 


2.87 


0.594 


0.509 


110.69 


15 


6.78 


0.580 


0.499 


106.17 


20 


12.65 


0.560 


0.483 


99.18 


25 


20.21 


0.534 


0.462 


90.12 


30 


29.88 


0.502 


0.434 


78.65 



(16) 



where % is the radius of the droplet (supposed to be spherical), 
and cr is the surface tension for the surface separating the hadron 
from the Q*-phase. The energy baiiier has a maximum at the 
critical radius 'Rc - '2-o-/[nQ.{fiH - f^Q')]- We neglected the term 
associated with the curvature energy, and also the terms con- 
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Fig. 3. Phase equilibrium curve for the hadron to Q*-matter 
phase transition. The continuous curve is relative to neutrino- 
free matter, the dashed curve to matter with trapped neutrinos. 
EOS: GMl plus MIT Bag model with B = 85 MeV/fm\ 
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B [MeV fm"^] 

Fig. 4. Phase equilibrium pressure Pq between the hadron phase 
and the Q* phase at T = 0, as a function of the bag constant 
B. The continuous curve refers to neutrino-free matter, and the 
dashed curve refers to neutrino-trapped matter 



nected with the electrostatic energy, since they are known to only 
introduce small corrections (lida & Sato 1998 Bombaci et al. 
I2004I I. The value of the surface tension cr for the interface sep- 
arating the quark and hadron phase is poorly known, and typi- 
cally values used in the literature range within 10 -50 MeV fm"^ 
(Heiselberg et al. 1993; lida & Sato 1998 ). We assume cr to be 
temperature independent and we take cr = 30 MeV fm"^. 

The quantum nucleation time can be straightforwardly 
evaluated within a semi-classical approach (Lifshitz & Kagan 
[1972 ; lida & Sato |1997i [1998). First one computes the ground 
state energy Eq in the WKB approximation and the oscillation 
frequency vo of the drop in the potential well UCR, T). Then, the 
probability of tunneling is given by 



Po = exp 



A(Eo) 



(17) 



where A(E) is the action under the potential barrier, which in a 
relativistic framework reads as (lida & Sato [1998) 



A(E) 



m{'R)c^ +E- UCRmUCR) - E] 



(18) 



with K-t the classical turning points and mCR) the droplet effec- 
tive mass. The quantum nucleation time is then equal to 



Tg = (VqPoNc) 



(19) 



with Nc ~ 10'^^ the number of nucleation centers expected in the 
innermost part (r < R„,ic ~ 100 m) of the hadronic star, where the 
pressure and temperature can be considered constant and equal 
to their central values. 

The thermal nucleation rate can be written (Langer & Turski 
[T973l )as 



— Qoexp(-f/(5R„r)/r) 
2n 



(20) 



where k is the so-called dynamical prefactor, which is related 
to the growth rate of the drop radius near the critical radius 
CRc)', i^o is the so-called statistical prefactor, which measures 
the phase-space volume of the saddle-point region around "Re', 
and UCRc, T) is the activation energy, i.e. the change in the free 
energy of the system required to activate the formation of a criti- 
cal size droplet. The Langer theory (Langer[l968] 119691 Langer 
& Turski 119731 Turski & Langer 1980 ) of homogeneous nu- 
cleation has been extended in Csernai & Kapusta (119921 ) and 
Venugopalan & Vischer (1994) to the case of first-order phase 
transitions occurring in relativistic systems, as in the case of the 
quark deconfinement transition. The statistical prefactor can be 
written (Csernai & Kapusta il992i) as 



fio 



(21) 



where is the quark correlation length, which gives a measure 
of the thickness of the interface layer between the two phases 
(the droplet "surface thickness"). In the present calculation, we 
take = 0.7 fm according to the estimate given in Csernai & 
Kapusta (1992) and Heiselberg ( 1 19951) . For the dynamical pref- 
actor, we have used a general expression, which has been derived 
by Venugopalan & Vischer ( 119941 ) 



2cr 



(22) 



where Aw = wq, - wh is the difference between the enthalpy 
density of the two phases, A the thermal conductivity, and tj and 
^ are the shear and bulk viscosities respectively of hadronic mat- 
ter The nucleation prefactor used in the present work differs sig- 
nificantly from the one used in previous works (Horvath et al. 
1992; Horvath [19941 Olesen & Madsen [l994b . where based on 
dimensional grounds, the prefactor was taken to be equal to . 

There are not many calculations of the transport properties of 
dense hadronic matter With a few exceptions (see e.g. van Dalen 
& Dieperink 2004; Chattaerjee & Bandyopadhyay 2006) most 
of them are relative to nuclear or pure neutron matter (Flowers 
&Itho [19791 D anielewicz [19841 Sedrakian et al. .l994. Benhar 
& Valli [20071 Shtemin & Yakovlev g008 , Benhar et al. [20T0l 
Zhang et al. '2010). These quantities have been calculated by 
Danielewicz 1984 in the case of nuclear matter According to 
the results of Danielewicz 1984 the dominant contribution to 
the prefactor k comes from the shear viscosity rj. Therefore, we 
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Fig. 5. Energy barrier for a virtual drop of the Q*-phase in 
stable neutrino-free hadronic matter as a function of the droplet 
radius and for different temperatures for a fixed pressure P - 
57 MeV/fml The bag constant takes the value B = 85 MeV/fm^ 

take A and ( equal to zero, and we use the following relation for 
the shear viscosity (Danielewicz 1984J : 

^ 7.6 xlO^S nH\2 MeV 
(r/MeV)2lno/ fms ' 

with «o = 0.16 fm"^ the saturation density of normal nuclear 
matter. 

The thermal nucleation time t^, relative to the innermost 
stellar region (V„uc = (47r/3)7?^„^) where almost constant pres- 
sure and temperature occur, can thus be written as 

T,/, = (y„„c/)"'. (24) 

In Fig. 5, we represent the energy barrier for a virtual drop 
of the Q*-phase (with B = 85 MeV/fm^) in the neutrino-free 
hadronic phase as a function of the droplet radius and for dif- 
ferent temperatures at a fixed pressure f = 57 MeV/fm^. As 
expected, from the results plotted in Fig. 1, the energy barrier 
UCR, T) and the droplet critical radius "Re decrease as the matter 
temperature is increased. This effect favors the Q*-phase forma- 
tion and, in particular, increases (decreases) the quantum nucle- 
ation rate (nucleation time r^) with respect to the corresponding 
quantities calculated at T = 0. 

In Fig. 6 we plot the quantum and thermal nucleation times 
of the Q*-phase (with B = 85 MeV/fm^) in jS-stable neutrino-free 
hadronic matter as a function of temperature and at a fixed pres- 
sure P - 51 MeV/fm^. As expected, we find a crossover tem- 
perature Tco above which thermal nucleation is dominant with 
respect to the quantum nucleation mechanism. For the case re- 
ported in Fig. 6, we have Tc, - 7.05 MeV and the corresponding 
nucleation time is log[o(T/s) = 54.4. The crossover temperature 
for different values of the pressure of y6-stable hadronic matter 
is reported in Table 4 (second column), together with the nucle- 
ation time calculated at T = Tco (third column). 

5. Evolution of proto-hadronic stars 

Keeping in mind the physical conditions in the interior of a 
proto-hadronic star (Burrow & Lattimer 1986? Prakash et al. 
[1997) (see Sect. 1 of the present paper) to establish whether 
this star will survive the early stages of its evolution with- 
out "decaying" to a quark star, one has to compare the quark 




T [MeV] 

Fig. 6. (Color online.) Thermal (r,/,) and quantum (r^) nucle- 
ation time of the Q*-phase (with B = 85 MeV/fm^) in yS-stable 
neutrino-free hadronic matter as a function of temperature at 
fixed pressure P = 51 MeV/fm^. The crossover temperature is 
Tco - 7.05 MeV. The limiting conversion temperature for the 
proto-hadronic star is, in this case, = 10.3 MeV, obtained from 
the intersection of the thermal nucleation time curve (continuous 
line) and the dot-dashed line representing log[Q(r/s) - 3. 



p 


T 

CO 


logio('^/s) 


53.98 


5.0 


233.6 


55.48 


6.0 


121.3 


56.94 


7.0 


56.6 


58.42 


8.0 


16.0 


58.85 


8.3 


3.0 



Table 4. Crossover temperature Tco (in MeV), for different fixed 
values of the pressure P (in Mev/fm^*) of neutrino-free hadronic 
matter The third column reports the logarithm of the nucleation 
time (in seconds) calculated at the crossover temperature. The 
value 8.3 MeV defines the value of the limiting conversion tem- 
perature for a star with a central pressure P - 58.85 MeV 
fm-l 



matter nucleation time t - mm{Tq, t,i,) with the cooling time 
tcooi ~ a few 10^ s. If T >> tcooi then quark matter nucleation 
is not likely to occur in the newly formed star, and this star will 
evolve into a cold deleptonized configuration. We thus introduce 
the concept of limiting conversion temperature for the proto- 
hadronic star and define it as the value of the stellar central tem- 
perature Tc for which the Q*-matter nucleation time is equal to 
10^* s. The limiting conversion temperature will clearly depend 
on the value of the stellar central pressure (and thus on the value 
of the stellar mass). 

The limiting conversion temperature is plotted in Fig. 7 as 
a function of the stellar central pressure P. A proto-hadronic star 
with a central temperature Tc > 0(f) will likely nucleate a Q*- 
matter drop during the early stages of its evolution, will finally 
evolve to a cold and deleptonized quark star, or will collapse to 
a black hole (depending on the particular model adopted for the 
matter EOS). 

For an isoentropic stellar core (Burrow & Lattimer 119861 
Prakash et al. 1 19971 1. the central temperature of the proto- 
hadronic star is given for the present EOS model by the lines la- 
beled by Ts in Fig. 7, relative to the case S - I ks and § -2 kg- 
The intersection point (Ps , 0s ) between the two curves 0(P) 
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35 40 45 
P [MeV fm " 

Fig. 7. (Color online.) The limiting conversion temperature © 
for a newborn hadronic star as a function of the central stellar 
pressure. The lines labeled Ts represent the stellar matter tem- 
perature as a function of pressure at fixed entropies per baryon 
§ Iks = 1 (dashed line) and 2 (solid line). Results for neutrino- 
free matter. 




105 120 135 150 165 
Pressure [MeV fm ] 

Fig. 8. (Color online.) The limiting conversion temperature for 
a newborn hadronic star as a function of the central stellar pres- 
sure for different values of the bag constant B. The lines labeled 
Ts represent the stellar matter temperature as a function of pres- 
sure at fixed entropies per baryon S Iks - 1, 1.5, 2. Results for 
neutrino-free matter. 



and Ts (P) thus gives the central pressure and temperature of the 
configuration that we denote as the critical mass configuration 
of the proto-hadronic stellar sequence. The value of the gravita- 
tional critical mass M„. = M{Ps , ©s) and baryonic critical mass 
MB,cr are reported in Table 5 for three different choices of the 
entropy per baryon, § /kg -Q (corresponding to a cold hadronic 
star)[[] , 1 and 2 in the case of a bag constant B = 85 Me V/fm^* and 



' I n Ref (Berezhiani et al. [20021 [20031 Bombaci et al. [20041 Drago 
et al. 120041 Lugones & Bombaci 2005; Bombaci et al. 120071 Bombaci 
et al. l2008t . the critical mass for cold (T = 0) metastable hadronic stars 
has been defined as the value of the gravitational mass for which the 
quantum nucleation time is equal to one year: M„-(r = 0) = M(t, = 
1 yr). It is worth recalling that the nucleation time is an extremely 
steep function of the hadronic star mass (Berezhiani et al. [20021 [20031 
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Fig. 9. (Color online.) Effects of neutrino trapping on the limit- 
ing conversion temperature © for a newborn hadronic star as a 
function of the central stellar pressure. 
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Fig. 10. (Color online.) Evolution of a proto-hadronic star in the 
gravitational-baryonic mass plane for B = 85 Me V/fm^ and cr = 
30 MeV/fm^. The upper (red) line represents the stellar equilib- 
rium sequence for neutrino-free proto-hadronic stars (PHS ) with 
S - Iks- The middle (blue) line represents the cold the HS se- 
quence. The asterisk and the full circle on these lines represent 
the stellar configuration with nucleation time t - oo and the ciit- 
ical mass configuration M„ , respectively. The lower (green) line 
represent the cold QS sequence. Assuming Mb - const, the evo- 
lution of a PHS in this plane occurs along a vertical line. Stellar 
masses are in units of the solar mass, M„,„ - 1.989 x lO^'-' e. 



for cr = 30 MeV/fm^. In the same table, we also report the value 
of the gravitational mass M of the cold hadronic star with bary- 
onic mass equal to Ms.cr- This configuration is stable (t = oo) 
with respect to Q*-matter nucleation in the case S Iks -2, and it 



Bombaci et al.'2004'), therefore the exact value of chosen in the def- 
inition of Mc,.{T = 0) is not crucial (one must take a "reasonable low" 
value of Tg, much shorter than the age of young pulsars as the Crab 
pulsar). We have verified that changing from I yr to 10^ s modi- 
fies M„(T = 0) by ~ 0.02%. On the other hand, the nucleation time 
T = min(T^,Trt) entering in the definition of the critical mass of proto- 
hadronic stars M„{S) must be comparable to the proto-hadronic star 
cooling time tcooi- 
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Fig. 11. (Color online.) Same as Fig. 10, but with the additional 
stellar sequence (upper (black) line) for proto-hadronic stars 
with trapped neutrinos (vPHSs) with a lepton fraction = 0.4 
and § -IkB- 




0.8 1 1.2 1.4 1.6 1.8 2 2.2 



Fig. 12. (Color onUne.) Same as Fig. 10, but with a bag constant 
B = 70.9 MeV/fml 



SIkB 


M„ 




M 


0.0 


1.573 


1.752 


1.573 


1.0 


1.494 


1.643 


1.485 


2.0 


1.390 


1.492 


1.361 



Table 5. Gravitational (M^r) and baryonic (Mb.ci) critical mass 
(see text for more details) for proto-hadronic stars at differ- 
ent entropy per baryon S Iks- M denotes the gravitational mass 
of the cold hadronic configuration with the same stellar bary- 
onic mass {MB,cr)- Stellar masses are in units of the solar mass, 
M,u„ = 1.989 X 10^^ g. Results are relative to a bag constant B 
= 85 MeV/fm^ 



is essentially stable (having a nucleation time enormously longer 
than the age of the universe) in the case S Iks - 1. 

The Umiting conversion temperature © for a newborn 
hadronic star is plotted in Fig. 8 for different values of the bag 
constant. The increase in B produces a growth of the region of 
the P-T plane where the proto-hadronic star could survive Q* 
nucleation and thus evolve to a cold hadronic star. 



To explore the role of neutrino trapping on the limiting con- 
version temperature and on the critical mass of the proto- 
hadronic star, we plot in Fig. 9 the results of our calculations for 
© and for the stellar matter temperature Ts at fixed entropies per 
baryon {§ /kg = 1, and 2) in the case of neutiino-free matter and 
neutrino-trapped matter. Results in Fig. 9 are relative to the GMl 
EOS for the hadronic phase and the bag model EOS with B = 85 
Mev/fm^ for the quark phase. As we see, neutrino trapping has a 
conspicuous effect on 0, and ultimately strongly suppresses the 
nucleation of Q*-matter in the proto-hadronic star. There is also 
a relatively small effect of neutrino trapping on Ts . 

The evolution of a proto-hadronic star (PHS) within our 
scenario is delineated in Fig. 10, where we plot the appropri- 
ate stellar equilibrium sequences in the gravitational-baryonic 
mass plane obained from the GMl EOS for the hadronic phase 
and the bag model EOS with B = 85 MeV/fm^ for the quark 
phase. In particular, we plot the PHS sequence, i.e. isoentropic 
HSs (5=2 ks) and neutrino-free matter (upper line), and the 
cold HS sequence (middle line). The asterisk and the full cir- 
cle on these lines identify respectively the stellar configuration 
with T = oo and the critical mass configuration. We denote as 
M™J = MB,cr(S - Iks) the baryonic critical mass for the PHS 
sequence and as M^^^. = Mb.cAS - 0) the baryonic critical mass 
for the cold hadronic star sequence. Finally, the lower line repre- 
sents the cold QS sequence. We assume (Bombaci & Datta 2000ll 
Mb - const during these stages of the stellar evolution ) Q. Thus 
according to the results in Fig. 10, proto-hadronic stars with a 
baryonic mass Mb < M^"J (- 1-492 Mj„„ within the present 
EOS parametrzation) will survive Q*-matter early nucleation 
(i.e. nucleation within the cooling time tcooi ~ a few 10^ s) and 
in the end will form stable (r = oo) cold hadronic stars. Proto- 
hadronic stars with M^^^ < Mb < A/g']',,^^. (the maximum bary- 
onic mass of the cold QS sequence, 1.813 M„,„ for the present 
EOS) will experience early nucleation of a Q*-matter drop and 
will ultimately form a cold deleptonized quark star. The last pos- 
sibility is for PHSs having Mb > '■^^^ ^^^^ early 
nucleation of a Q*-matter drop will trigger a stellar conversion 
process to a cold QS configuation with Mb > A/f thus these 
PHSs will finally form black holes. 

The outcomes of this scenario are not altered by neutrino 
trapping effects in hot y6-stable hadronic matter as illustrated 
in Fig. 1 1 where (in addition to the stellar sequence curves of 
Fig. 10) we plot the stellar sequence for proto-hadronic stars 
with trapped neutrinos (vPHSs) with a lepton fraction Yi - 0.4 
and S - 2 kB- For the vPHS stellar sequence in Fig. 11 we 
find a critical baryonic mass M^™'^ - 1.96 Msm- Thus vPHSs 
with Mb > Mg™^ after neutrino escape and cooling will finally 
evolve to black holes. The fate of a vPHS with Mb < M^g"^ 
is the same as the corresponding neutrino-free PHS with equal 
baryonic mass. 

Finally, in Fig. 12 we plot the PHS, cold HS, and cold QS 
sequences in the gravitational-baryonic mass plane for the case 
of a bag constant B = 70.9 MeV/fml (All the other EOS pa- 
rameters are the same as for the results in Fig. 10.) Apart from 
changes in the numerical values of the maximum masses of the 
various stellar sequences and of the values of the stellar criti- 
cal masses, the evolutionary scenario of a proto-hadronic star is 



- Sizeable mass accretion on the proto-neutron star occurs within a 
time of ~ 0.5 s after core bounce. (Burrow &. Lattimer ll986l Prakash 
et al. I1997I I. During the subsequent stages, the star thus evolves with 
Mb - const. 
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qualitatively similar to the one previously discussed for the re- 
sults in Fig. 10. 

An interesting situation comes up when a PHS with Mb < 
M^^^ is formed in a binary stellar system. In this case, as we 
have just seen, a stable (case of Fig. 10) or a metastable (case 
of Fig. 12) cold HS is formed. Long-term accretion from the 
companion star will start to populate the portion of the cold HS 
sequence with M^^^ < Mb < M^^^^. If the star accretes enough 
matter (~ 0.1 - 0.3 Ms„„, see Figs. 10 and 12), it could possibly 
reach the critical mass M^^^.. At this point a quantum nucleation 
process will trigger the conversion of this critical mass cold HS 
to a QS (Berezhiani et al. l2002l l2003l Bombaci et al. (2004l l. 
The time delay between the supernova explosion forming the 
PHS, and the second "explosion" (neutrino burst or/and GRB) 
forming the QS will actually depend on the mass accretion rate, 
which modifies the quark matter nucleation time via an explicit 
time dependence of the stellar central pressure. 

Within the present scenario, there is thus a stellar mass range 
Mg"^ < Mb < Mgl^ where it is possible to find two differ- 
ent types of "neutron stars" in the Universe: pure hadronic com- 
pact stars with large radii in the range 12-20 km, and quark 
stars with small radii in the range of 6 - 9 km (Bombaci et 
al. [2004; Drago & Lavagno 120 1 Oi l. Accurate measurements of 
both the mass and radius of a few individual "neutron stars" 
(Bhattacharyya 2010; Steiner et al. 2010) could shed light on 
the validity of this presumptive scenario. 



6. Limiting mass of compact stars 

The possibility of having metastable hadronic stars, together 
with the feasible existence of two distinct families of compact 
stars (pure hadronic stars and quark stars), demands an extension 
of the concept of maximum mass of a "neutron star" with re- 
spect to the classical one introduced by Oppenheimer & VoUcoff 
(1939). Since metastable HS with a "short" mean-life time are 
very unlikely to be observed, an extended concept of maximum 
mass has to be introduced in view of the comparison with the 
values of the mass of compact stars deduced from direct astro- 
physical observation. With this operational definition in mind, 
Bombaci et al. (2004) define the limiting mass {Mum) of a com- 
pact star in the case of cold stellar configurations. This concept 
of Umiting mass can be straightforwardly extended to the present 
case, i.e. taking the effects of proto-neutron stars evolution into 
account. For the two cases reported in Figs. 10 and 12, the lim- 
iting mass is the highest between the gravitational critical mass 
M^'^ of the cold HS sequence and the gravitational maximum 
mass M^^ax of the cold QS sequence. Thus M/,„, - M^^ (case in 
Fig. 10) and Mum - M^,l^ (case in Fig. 12). 

The very recent measurement (Demorest et al. 12010^ of the 
Shapiro delay in the binary millisecond pulsar J1614-2230 has 
helped to obtain the mass of the associated neutron star. The 
calculated mass is M = (1.97 ± O.O4)M0 (Demorest et al. 120 1 Oi l, 
making PSR J1614-2230 the most massive neutron star known 
to date. 

As is well known (see e.g. Lattimer & Prakash 120011 
Bombaci 2007), neutron star mass measurements give one of the 
most stringent tests of the composition and EOS of strong in- 
teracting matter at very high densities. Here, we very briefly ad- 
dress some possible implications of the high mass of PSR J1614- 
2230 in connection with the feasibility of the scenario discussed 
in present work. To this aim we consider two almost diametrical 
situations for the stiffness of the EOS. 



The first possibility is for the case of a stiff hadronic EOS 
and a soft quark matter EOS. This is the case, for example, for 
the GMl model with jCq- - Xp - 0.8, x„ = 0.913 for hyperonic 
matter, and the MIT bag model EOS (Farhi & Jaffe 1984)) with 
B - 100-150 MeV/fm^ for quark matter For these EOS mod- 
els, the limiting mass is given by the gravitational critical mass 
of the HS sequence, and it is in the range M/„„ = 1.9-2.1M0 
(Bombaci et al. l2008l l. Thus, in this case, PSR J1614-2230 can 
be interpreted as a pure hadronic star (hyperon star), and quark 
matter nucleation, in HS with M > M^'^, will produce stellar 
configurations which will form black holes. 

The second possibility is for a soft hadronic EOS and a 
stiff quark matter EOS. This is the case, for example, of the 
microscopic Brueckner-Hartree-Fock EOS for hyperonic mat- 
ter (Baldo et al. l2000l Vidana et al. '2000', Schulze et al. l2006b 
and the recent perturbative calculations of the quark matter EOS 
by (Kurkela et al. 120 1 Oi l which include quark interaction ef- 
fects up to the second order in the QCD coupling Qs (see also 
Fraga et al. 120011 Alford et al. l2005T l. For these EOS models, 
the limiting mass is likely given by the gravitational maximum 
mass of the cold QS sequence, which according to the results of 
(Kurkela et al. l20T0] l is M,^f , ~ 2Mq (hybrid stars), M^L = 2.0- 
2.7Mq (strange stars). Thus, for this second case, PSR J1614- 
2230 could be interpreted as a quark star (Ozel et al. 120 1 Oi l. 

7. Conclusions 

In summary, in this work we have studied the quark deconfine- 
ment phase transition in hot y6-stable hadronic matter and ex- 
plored some of its consequences for the physics of neutron stars 
at birth. We calculated and compared the nucleation time due 
to thermal and quantum nucleation mechanisms, and computed 
the crossover temperature above which thermal nucleation dom- 
inates the finite temperature quantum nucleation mechanism. In 
addition, we introduced the new concept of limiting conversion 
temperature for proto-hadronic stars and extended the con- 
cept of critical mass (Berezhiani et al. 120021 120031 Bombaci et 
al. l2004l l to the case of finite temperature hadronic stars. 

Our main finding is that proto-hadronic stars with a gravita- 
tional mass lower than the critical mass M„ could survive the 
early stages of their evolution without decaying to a quark star. 
This outcome contrasts with the predictions of the earlier stud- 
ies (Horva flietal . [T992l Horv adilT994l Olesen & Madsen [T994l 
Heiselberg |I995l Harko et al. '2004) where it was inferred that 
all the pure hadronic compact stars, with a central temperature 
above 2-3 MeV, are converted to quark stars within the first 
seconds after their birth. 

However, the prompt formation of a critical size drop of 
quark matter could take place when M > Mcr- These proto- 
hadronic stars evolve to cold and deleptonized quark stars or 
collapse to a black holes. 

Finally, if quark matter nucleation occurs during the post- 
bounce stage of core-collapse supernova, then the quark decon- 
finement phase transition could trigger a delayed supernova ex- 
plosion characterized by a peculiar neutrino signal (Sagert et al. 
120091 Mintz et al. ilOTOl Nakazato et al. l2008l Dasgupta et al. 
I20T01) . 
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